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Dynamics of p53 is known to play important roles in the regulation of cell fate decisions in response 
to various stresses, and PDCD5 functions as a co-activator of p53 to modulate the p53 dynamics. 

In the present paper, we investigate how p53 dynamics are modulated by PDCD5 during the DNA 
damage response using methods of bifurcation analysis and potential landscape. Our results reveal 
that p53 activities can display rich dynamics under different PDCD5 levels, including monostability, 
bistability with two stable steady states, oscillations, and co-existence of a stable steady state and 
an oscillatory state. Physical properties of the p53 oscillations are further shown by the potential 
landscape, in which the potential force attracts the system state to the limit cycle attractor, and 
the curl flux force drives the coherent oscillation along the cyclic. We also investigate the effect of 
PDCD5 efficiency on inducing the p53 oscillations. We show that Hopf bifurcation is induced by 
increasing the PDCD5 efficiency, and the system dynamics show clear transition features in both 
barrier height and energy dissipation when the efficiency is close to the bifurcation point. This 
study provides a global picture of how PDCD5 regulates p53 dynamics via the interaction with the 
p53-Mdm2 oscillator and can be helpful in understanding the complicate p53 dynamics in a more 
complete p53 pathway. 


I. INTRODUCTION 

The tumor suppressor p53 plays a central role in cellu¬ 
lar responses to various stress, including oxidative stress, 
hypoxia, telomere erosion and DNA damage [lj [2j. In 
unstressed cells, p53 is kept at low level via its neg¬ 
ative regulator Mdm2 [3]. Under stressed conditions, 
such as DNA damage, p53 is stabilized and activated 
to induce the express of downstream genes, including 
p21/WAFl/CIPl and GADD45 that are involved in cell 
cycle arrest, and PUMA, Bax and PIG3 that can in¬ 
duce apoptosis ms]. The cell fate decisions after DNA 
damage are closely related to the p53 dynamics that is 
regulated by p53-Mdm2 interactions 0 m. Oscillations 
of p53 level have been observed upon IR induced DNA 
damage at the population level in several human cell lines 
and transgenic mice 0EHIO]. More interestingly, pulses 
of p53 level were revealed in individual MCF7 cells, and 
it was suggested that the cell fate is governed by the 
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number of p53 pulses, i.e., few pulses promote cell sur¬ 
vival, whereas sustained pulses induce apoptosis HU eng. 
Fine control of p53 dynamics is crucial for proper cellular 
response. Mutations and deregulation of p53 expression 
have been found to associate with various cancer types 

mm- 

Programmed Cell Death 5 (PDCD5), formerly referred 
to as TFAR19 (TF-1 cell apoptosis-related gene 19), is 
known to promote apoptosis in different cell types in re¬ 
sponse to various stimuli m- Decreased expression of 
PDCD5 has been detected in various human tumors m 
! l8L and restoration of PDCD5 with recombinant protein 
or an adenovirus expression vector can significantly sensi¬ 
tive different cancers to chemotherapies HUEE9]. PDCD5 
is rapidly upregulated after DNA damage, interacts with 
the apoptosis pathway, and translocates from the cyto¬ 
plasm to nucleus during the early stages of apoptosis 
[20l [22] . Recently, novel evidence indicates that PDCD5 
is a p53 regulator during gene expression and cell cy¬ 
cle [23] . It was shown that PDCD5 interacts with p53 
by inhibiting the Mdm2-mediated ubiquitination and ac¬ 
celerating the Mdm2 degradation. Hence, upon DNA 
damage, PDCD5 can function as a co-activator of p53 to 
regulate cell cycle arrest and apoptosis. 

A series of computational models have been con- 
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structed to investigate the mechanism of p53-mediated 
cell-fate decision ia 0 in mm no. In these models 
the p53-Mdm2 oscillation was considered to be crucial 
for cell-fate decision after DNA damage. After DNA 
damage, such as double strand breaks (DSBs), active 
ATM monomer (ATM*) become dominant. In the nu¬ 
cleus, ATM* active p53 in two ways: ATM* promotes 
the phosphorylation of p53 at Ser-15 [26. and acceler¬ 
ates the degradation of Mdm2 in nucleus (Mdm2 nuc ) [27j . 
In the cytoplasm, active p53 induces the production of 
Mdm2 cyt , which in turn promotes the translation of p53 
mRNA to form a positive feedback. The p53-Mdm2 oscil¬ 
lation is a consequence of the two coupled feedback loops: 
the negative-feedback between p53 and Mdm2 nuc , and 
the positive-feedback between p53 and Mdm2 cyt . Re¬ 
cently, mathematical models of how PDCD5 interacts 
with the DNA damage response pathway to regulate cell 
fate decisions have been developed [28, .29 . In our pre¬ 
vious study [29] . a model for the effect of PDCD5 to the 
p53 pathway has been established through a nonlinear 
dynamics model that is developed based on the mod¬ 
ule of p53-Mdm2 oscillator in mm and experimental 
findings in [23]. It was shown that the p53 activity can 
display different dynamics after DNA damage depend¬ 
ing on the PDCD5 level. The p53 protein shows low 
activity in case of PDCD5 deletion, sustain intermediate 
level for medial PDCD5 expression, and pulses when the 
PDCD5 level is upregulated [29]. Nevertheless, little is 
known about the global p53 dynamics upon PDCD5 in¬ 
teractions with changes in the expression levels of p53 
and PDCD5, which is often seen in tumors. 

Here, we investigate, based on the mathematical model 
proposed in m, how the dynamics of p53 activity after 
DNA damage depend on changes in the levels of p53 pro¬ 
duction and PDCD5. A global bifurcation analysis shows 
that p53 activity can display various dynamics depend¬ 
ing on the p53 production and PDCD5 levels, including 
monostability, bistability with two stable steady states, 
oscillations, and co-existence of a stable steady state and 
an oscillatory state. These dynamics are further investi¬ 
gated through the method of potential landscape. The 
stability of the oscillation states are characterized by the 
potential force and the probability flux. We also discuss 
the effect of PDCD5 efficiency on p53 dynamics. PDCD5 
efficiency can induce p53 oscillation by Hopf bifurcation, 
and the transition of p53 dynamics near the Hopf bi¬ 
furcation are charaterized by barrier height and energy 
dissipation. 


II. MODEL AND METHOD 

A. PDCD5-mediated p53-Mdm2 oscillator 

Figure, [l] illustrates the model of p53-Mdm2 oscillator 
with PDCD5 regulation studied in this paper. Here we 
summarize the model equations and refer [29] for details. 
The model equations are composed of three compo¬ 



FIG. 1. Illustration of the PDCD5-mediated p53-Mdm2 path¬ 
way. Active p53 promote the production of Mdm2 cyt , which 
in turn promote the translation of p53 mRNA. In the nu¬ 
cleus, active p53 is degraded by binding to Mdm2 nU c, and the 
interaction is disrupted by PDCD5. Both active ATM and 
PDCD5 are able to accelerate the degradation of Mdm2 nU c- 
Mdm2 in the nucleus and cytoplasm can shuttle between the 
two compartments. Refer to the text and [29] for details. 


nents: active p53 in the nucleus [p53], Mdm2 in the 
nucleus [Mdm2 nuc ] and in cytoplasm [Mdm2 cyt ]. After 
DNA damage, ATM is activated and promotes the degra¬ 
dation of Mdm2 nuc so that the p53 level is unregulated. 
Active p53 promotes the production of Mdm2 cyt , which 
promotes the translation of p53 mRNA to produce p53 
and hence form a positive-feedback loop. In the nucleus, 
active p53 is degraded slowly at a rate d P 53 by weakly 
binding to Mdm2 nuc , and the interaction is disrupted 
by PDCD5. Moreover, PDCD5 can also accelerate the 
degradation of Mdm2 nuc . Mdm2 in the nucleus and cy¬ 
toplasm can shuttle between the two compartments at 
rates k- in and k Qut , respectively. These interactions give 
the following differential equations: 

v p53 ([Mdm2 cyt ]) - d p53 ([Mdm2 nuc ]) [p53] 

Fi([p53], [Mdm2 cyt ], [Mdm2 nuc ]), (1) 

^Mdm 2 ([p53]) - k in [Mdm2 cyt ] 

k ou t [Mdm2 nuc ] dMdm 2 [Mdm2 cy t] 
F 2 ([p53], [Mdm2 cyt ], [Mdm2 nuc ]), (2) 

k in [Mdm2 cyt ] - k out [Mdm2 nuc ] 

- fit) d M dm 2 [Mdm2 nuc] 

C 3 ([p53], [Mdm2 cyt ], [Mdm2 nuc ]). (3) 


Here /q n , & 0 ut> ^Mdm 2 are constants, and other rate func- 


d[p53] 

dt 

d[Mdm2 cyt ] 

dt 


d[Mdm2 nuc ] 

dt 
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TABLE I. Typical parameter values (The unit for time is min, 
and arbitrary unit for the concentration) (Ref. 29]). 


Si 

4 

Pi 

0.991 

Ki 

0.057 

Vp53 

0.85 

S2 

4 

p2 

0.9873 

dp53 

0.4 

m\ 

4 

ri 

0.8 

on 

3.3 

k 2 

0.09 

S3 

4 

P3 

0.989 

K 3 

4.433 

t’Mdm2 

0.135 

kin 

0.14 

kout 

0.01 

^Mdm2 

0.036 

f 

2.7 

S4 

4 

P4 

0.2 

Ka 

0.41 

S5 

4 

P5 

0.5 


1.58 






tions are given by 


i’ P 53 ([Mdm 2 cyt ]) = v p53 (1 - Pi ) + Pi 


[Mdm2 


cytj 


K? + [Mdm2 cyt ] 


^p 53 ([Mdm2 nuc ]) — <ip 53 (1 P 2 ) 


+ p2 


[Mdm2 n 


^Mdm2 ([p53]) = ^Mdm2 (1 “ P?>) + P3 

fit) = f\(l-p4- P5) + P4 


K 2 {P{t )) s2 + [Mdm 2 nuc ] S2 J 
[p53p 


K, 3 s3 + [ P 53p 

p(ty 4 

K[ 4 + P(t)«‘ 


+ P5 


A(t) 


s 5 


K s y + A(ty* 


where 


K 2 (P) = K 2 (( 


( 4 ) 


Here A(£) and P(t) are time dependent functions for the 
levels of active ATM and PDCD5 in nucleus, respectively. 

In this study, we only consider the process of DNA re¬ 
pair after DNA damage that can persist for as long as 
10 hours before the onset of apoptosis [22 . Both active 
ATM and PDCD5 level are upregulated after DNA dam¬ 
age. Hence, we set A(t) = 5 to mimic this process in 
accordance with the simulation given by mm- Simi¬ 
larly, during DNA repair, we assume P(t) = Po with Po 
represents the PDCD5 expression level. 

We refer Table [I] for parameter values in this study. 


of the non-equilibrium interactions. Here we summarize 
the potential landscape and flux framework that were in¬ 
troduced by Wang et al. in [33H35] . 

We write the dynamical equation of the p53- 
Mdm2 oscillator as X = F(X), where X = 

([p53], [Mdm2 cyt ], [Mdm2 nuc ]) and F(X) represents the 
right hand size of Q- <©• The above equation can be ex¬ 
tended to include fluctuations by a probability approach: 
X = F(X) + £, where ( is the noise perturbation. The 
statistical nature of the noise is often assumed as Gaus¬ 
sian (large number theorem) and white (no memory): 
(CCOCCO) = 2D S(t — t') and (£(£)) = 0, where D is the 
correlation tensor (matrix) measuring the noise strength. 

The probability of system states P(X, t) evolves follow¬ 
ing the Fokker-Plank equation ^+V-J(X, t) = 0, where 
the flux vector J is defined as J(X, t) = F P — D • J^P. 
The flux J measures the speed of the flow in the concen¬ 
tration space. 

At the stationary state, ^ = 0, then V • J(X,£) = 0. 
There are two possibilities: one is J = 0, which implies 
detailed balance, and the other is J ^ 0 so that the 
detailed balance is broken and the system is at the non¬ 
equilibrium state. 

At the stationary state, the probability flux is defined 
by 

J ss = F P ss - D • -t^Pss (5) 

where P ss is the probability density at stationary state. 
Hence, 


F — D • T^rPss/ P ss + Jss/Pss 
= D ' 1 nP ss ) + Jss/Pss 

= --D--^U + 3 ss /P ss . ( 6 ) 

Here U = — In P ss is defined as the non-equilibrium po¬ 
tential. From (J 6 |, we have divided the force F into two 
parts: the potential force —D • -§^U and curl flux force 
J ss /Pss- At detail balance, the curl flux force is zero. For 
a non-equilibrium open system, both potential landscape 
and the associated flux are essential in characterizing the 
global stationary state properties and the system dynam¬ 
ics. 


B. Potential landscape and probability flux 

The ideas of potential landscape have been introduced 
for uncovering global principles in biology for protein 
dynamics [30] . interactions mm, and gene networks 
[33H35] . For a non-equilibrium open systems such as the 
p53-Mdm2 oscillator that exchanges energies and infor¬ 
mations with outside environments, the potential and 
non-equilibrium probability flux are able to reveal in¬ 
sights for the global robustness and physical mechanisms 


III. RESULTS 

A. Codimension-two bifurcation analysis 

To investigate how PDCD5 interacts with p53 to reg¬ 
ulate cell fate decision dynamics after DNA damage, we 
performed codimension-two bifurcation analysis with re¬ 
spect to the two parameters Po and ^ 53 . The bifurca¬ 
tion diagrams were computed with AUTO incorporated 
in XPPAUT [36]. The main bifurcation diagram is shown 
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at Fig. [2] and detailed below, with Fig. It the enlarge of 
the dashed square region in Fig. It- 

The parameter plane (^ 53 , Po) is divided into ten re¬ 
gions labelled I-X each with different dynamical profiles 
and marked with different colors in Fig. [ 2 ] (refer Fig. [2J) 
for regions VIII-X). Typical dynamics of each region is 
shown at Fig. [3] and described below: 

1. Region I corresponds to monostability with a single 
stable steady state of low p53 activity despite the 
expression level of PDCD5. 

2. In region II there is a stable steady state and two 
unstable steady states. Regions I and II are sep¬ 
arated by a curve /* of saddle-node bifurcation, 
across which an unstable node and a saddle appear 
due to fold bifurcation. 

3. Region III is on the right of region II. Crossing the 
curve snic from region II, the stable node and the 
saddle collide and disappear, accompanied by the 
emergence of a stable limit cycle that surrounds an 
unstable steady state. Meanwhile, parameters can 
also change from region I to III crossing the bound¬ 
ary suphi or subhi, respectively (Fig. [2p. The 
curve suphi represents supercritical Hopf bifurca¬ 
tion at which the stable steady state in region I 
becomes unstable and a stable limit cycle appears. 
The curve subhi represents subcritical Hopf bifur¬ 
cation. Before crossing the curve subhi, two limit 
cycles, stable and unstable ones, come out due to 
fold bifurcation of limit cycles that is represented 
by the curve Ipci (Fig. It)- Next, the unstable 
limit cycle disappear and the stable steady state 
lose its stability when at the curve subhi. 

4. Regions III and IV are separated by a subcritical 
Hopf bifurcation curve subh 2 by which an unstable 
limit cycle appears and locates between the stable 
limit cycle and the steady state, and the steady 
state that is unstable at region III becomes stable 
at region IV. 

5. Region V, similar to region I, has a stable steady 
state which corresponds to a state of high p53 con¬ 
centration. The two limit cycles at region IV col¬ 
lide and disappear when parameters cross the fold 
bifurcation curve lpc 2 that separates regions IV 
and V. 

6. Three steady states arise in region VI, two of them 
are unstable and appear while crossing the fold bi¬ 
furcation curve f 2 of equilibria between regions VI 
and V. 

7. Region VII gives bistability in p53 activity, in¬ 
cluding two stable steady states with low or high 
p53 concentrations, and an unstable steady state. 
There are two ways to reach this bistable region. 
From region VI and crossing the subcritical Hopf 
bifurcation subhi, an unstable limit cycle emerges 


and the unstable steady state with low p53 level 
becomes stable, but the unstable limit cycle im¬ 
mediately disappear due to homoclinic bifurcation 
(refer the curve hci in Fig. |2^i which is very close 
to subhi). From region / and crossing the bound¬ 
ary /{“, a pair of stable and unstable steady states 
appear due to fold bifurcation of equilibria. 

8. There are three regions ( VIII-X) around 
(f^p 535 To) = (1-06,0.41) which share the same 
phase diagram with a stable limit cycle and three 
steady states (Fig. [ 2 ) 3 ) . In region VIII the 
three steady states are unstable, two of them (a 
saddle and an unstable focus) come from the fold 
bifurcation (/+) that separates regions III and 
VIII. In region IX the steady state with high 
p53 level becomes stable and is surrounded by an 
unstable limit cycle. The change is originated from 
the subcritical Hopf bifurcation subh 2 from region 
VIII. Finally, the unstable limit cycle disappear 
in region X due to the homoclinic bifurcation hc 2 
between regions IX and X. We note that in regions 
IX and X, there are bistable states of a stable 
steady state and a stable limit cycle, similar to the 
region IV. 

In addition to the above ten regions, there are four 
codimension-two bifurcation points denoted by black 
dots in Fig. [ 2 ] two cusp points (CPi and CP 2 ), one 
Generalized Hopf bifurcation (GH) and one Bogdanov- 
Takens bifurcation (BT). The cusp point CPi locates 
at (^ P 53 ,Po) = (0.6464,0.508) where the fold bifurca¬ 
tion curve /* and the saddle-node homoclinic bifurcation 
curve snic meet tangentially. The other cusp point CP 2 
is given by (1.0220,0.4267) by the two fold bifurcation 
curves fi and , where two fold bifurcations coalesce and 
disappear. The point GH at (1.0071,0.4) corresponds to 
the meeting point of suphi and subhi from where a fold 
bifurcation curve of the limit cycle (Ipci) occurs. At the 
BT point (1.0885,0.3873) the fold bifurcation curve fi 
and subh 2 meet tangentially to give a homoclinic bifur¬ 
cation curve hc 2 . The BT point separates fi into two 
segments /{" and f+. Crossing /-f from left to right 
gives a stable equilibrium and a saddle, while crossing 
fi gives an unstable equilibrium and a saddle m- 

In summary, by manipulating the expression level of 
PDCD5 and the maximum production rate of p53, the 
system displays four types of stable dynamics: a single 
stable steady state (regions I, II, V, and VI), two stable 
steady states (region VII), a stable limit cycle (regions 
III and VIII), and coexistence of a stable steady state 
and a stable limit cycle (regions IV, IX, and X). The 
bifurcation diagram is zoomed out by codimension-one 
bifurcations with seven Pq values marked by a-g in Fig. 
12 ] and are detailed at the next section. 

From Fig. [2] during DNA damage and when the 
PDCD5 is upregulated (Pq > 0.3873), the system shows 
either low p53 activity, p53 oscillation, or sustained high 
p53 level with the increasing of p53 production rate. p53 
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FIG. 2. Codimension-Two bifurcation diagram, (a) The two-parameter bifurcation diagram with respect to u P 53 and Po- 
(b) Enlarged representation of the rectangle in (a). Codimension-one bifurcation curves are: suph -the supercritical Hopf 
bifurcation; subhi, subh 2 ~ the subscribe hopf bifurcation; p-the fold bifurcation of equilibria; lpc\, lpc 2 ~fo\d bifurcation of limit 
cycles; /ici, hc 2 -homoclinic bifurcation. These curves mainly divide the (u P 53 ,Po) plane into ten regions I-X . Codimension- 
two bifurcation points are: GH-generalized Hopf bifurcation, CP-the cusp bifurcation and BT-Bogdanov-Takens bifurcation. 
Values of Po denoted by a-g correspond to Po values in Fig. [5] respectively. 


oscillations have been known to be essential for cellu¬ 
lar response to DNA damage including DNA repair and 
apoptosis ® nn m [39] . Our analyses suggest that 
proper cell response to DNA damage with p53 oscillation 
is possible only when PDCD5 is highly expressed and the 
p53 production rate takes proper value. To further inves¬ 
tigate the oscillation dynamics, Fig. [4] shows the period 
and amplitude of the stable limit cycles corresponding to 
regions III and IV in Fig. [ 2 ] The results show that the 
oscillation periods decrease with the maximum p53 pro¬ 
duction rate (^ 53 ), while the amplitudes increase with 
F P 53 . This is consistent with the experimental observa¬ 
tions m • However both periods and amplitudes are in¬ 
sensitive with the PDCD5 level over a wide parameter 
range. These results indicates that upregulated PDCD5 
and a proper p53 expression level are essential for pro¬ 
ducing proper p53 oscillation that regulates the cellular 
response to DNA damage. 

When PDCD5 is downregulated (Po < 0.3873), the 
system has either a single steady state with low p53 activ¬ 
ity or bistability with either low or high p53 activity after 
DNA damage depending on the maximum p53 produc¬ 
tion rate. Experiments have shown that cells exposed to 
sustained p53 signaling frequently underwent senescence 

m • 

B. Codimension-one bifurcation analysis 

To get a clear insight into the codimension-two bifur¬ 
cation diagram in Fig. [2] we considered a codimension- 


one bifurcation of the concentration of p53 with respect 
to the p53 production rate v p ^. The bifurcation dia¬ 
gram is shown at Fig. [5] with panels a-g for given values 
Po = 0.3,0.39,0.41,0.42,0.48,0.51,0.8, respectively, as 
marked at Fig. [ 2 ] Diagrams for each Po value is detailed 
below. 

(a) . For low PDCD5 level (Po = 0.3), the equilibrium 

show a S-shaped bifurcation diagram with bistability 
(Fig. [5p. There are three branches of the equilib¬ 
rium depending on ^ 53 . The upper branch is com¬ 
posed of stable nodes, states at the middle branch 
are saddles, while the lower branch consists of stable 
nodes and foci which are separated near the subcrit- 
ical Hopf bifurcation point Hi ( v p 53 = 2.315). There 
is a fold bifurcation of equilibria Pi (Fp 53 = 1.143) 
where a stable node and a saddle appear, so that 
bistability occurs with F P 53 between Pi and Hi. 
With the increasing of ^ 53 , an unstable limit cycle 
appears from a homoclinic bifurcation point HCi at 
a saddle (Fp 53 = 2.239), and the unstable limit cy¬ 
cle disappears via the subcritical Hopf bifurcation 
Hi. Crossing Pi, stable focus at the lower branch 
change to unstable one. The saddle and the lower 
unstable focus collide and disappear via another fold 
bifurcation P 2 at ^53 = 3.305. 

(b) . When Po increases (Po = 0.39) (Fig. [5 ]d), the fold 

bifurcation Pi becomes larger than the subcritical 
Hopf bifurcation Hi so that the bistability of equilib¬ 
rium vanishes. With the increasing of Fp 53 , an unsta¬ 
ble limit cycle arises from the Hopf bifurcation Hi. 
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FIG. 3. Phase diagrams corresponding to the labeled regions I-X in Fig. [2] Red solid dots represent stable equilibria, black open 
circles unstable one; green solid and blue dash dotted lines represent stable limit cycles and unstable limit cycles respectively. 



0.6 0.8 1 0.6 0.8 1 

q>53 ^p53 


FIG. 4. Periods (a) and amplitudes (b) of oscillations corre¬ 
sponding to the regions of oscillation dynamics in Fig. [2] 


The limit cycle becomes stable at the fold bifurcation 
LPC \, and then vanishes at the second fold bifurca¬ 


tion LPC 2 . Different from the case Pq = 0.3, now 
the fold bifurcation F\ gives a saddle (the middle 
branch) and an unstable focus (the upper branch) 
(Fig. pb inset). The unstable focus on the upper 
branch becomes stable and an unstable limit cycle 
appears via the subcritical Hop bifurcation H 2 . The 
unstable limit cycle develops until it meet the sad¬ 
dle at the homoclinic bifurcation point HC 2 . The 
stable focus and stable limit cycle coexist with 
between H 2 and LPC 2 , which gives the region X in 
Fig.0 

(c) . When Pq increases to 0.41 (Fig. [5]:), the diagram 

is similar but the Hopf bifurcation Hi changes from 
subcritical to supercritical, so that stable limit cycles 
exist with n P 53 taken values over a wide range from 
Hi to LPC 2 . 

(d) . When P 0 = 0.42, Pi becomes too close to P 2 so that 

the unstable limit cycle arise from H 2 does not col- 
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lide with a saddle but ends up at the fold bifurcation 
point of limit cycles LPC 2 (Fig. m 

(e) . For a larger Po (Fig. §0 , the two fold bifurcation 

points Fi and F\ coalesce and disappear due to a 
codimension-two cusp point CP 2 in Fig. [2] 

(f) . With the further increasing of Po (Fig. [5]f), another 

two fold bifurcation points (F\* and P 2 ) arise from 
the Hopf bifurcation Hi due to the cusp point CP\ 
in Fig. [2] 

(g) . Further increasing Po to higher level (Po = 0.8), the 

supercritical Hopf bifurcation point Hi moves to the 
left and collides with P 2 * to give a fold-Hopf bifurca¬ 
tion point ( ZH), and then P 2 changes to a saddle- 
node homoclinic bifurcation point SNIC, from which 
the saddle and the node collide and disappear to gen¬ 
erate a stable limit cycle. 

The above bifurcation analyses present the global dy¬ 
namics of the whole system, which is further explored be¬ 
low through the underlying potential landscape [30] Up. 


C. Potential landscapes and global dynamics 

To explore the global dynamics from the potential per¬ 
spective, we projected the potential function to two in¬ 
dependent variables, the p53 concentration [p53], and 
the total Mdm2 concentration [Mdm2] = [Mdm2 cyt ] + 
[Mdm2 nuc ]. The potential landscapes for parameter val¬ 
ues taken from the 10 typical regions are shown at Fig. 
[ 6 ] 

The potential landscapes show that when there is a 
single stable steady state (I, II, V, VI), the potential is 
funnelled towards a global minimum which corresponds 
to the global stable steady state. 

In region VII, there are two stable steady states, and 
the potential has two local minimum, corresponding to 
high or low p53 levels, respectively. The landscape at low 
p53 state has wide attractive region and shallow slop, 
while at the high p53 state has small attractive region 
and deep slope. These suggest that a cell from randomly 
select initial state is more likely to response with low p53 
state because of the wider attractive region. 

In other regions with p53/Mdm2 oscillations (III, IV, 
VIII, IX, X ), the potential landscapes show an irregular 
and inhomogeneous closed ring valley that corresponds to 
the deterministic stable limit cycle trajectory. It is obvi¬ 
ous that the landscape is not uniformly distributed along 
the limit cycle path due to the inhomogeneous speed on 
the limit cycle mim. We also note that in the land¬ 
scapes for regions IV, IX and X, in addition to the closed 
valley, there is a deep funnel towards a local minimum. 
These local minimum show the coexistence of the stable 
steady state (marked by P and the arrows in Fig. [ 6 |. 

The regions with p53 oscillations are the most inter¬ 
ested because the oscillation dynamics are essential for 


cell fate decision in response to various stresses. To fur¬ 
ther analyze the potential landscape when the system 
display stable p53 oscillations, we examined the case at 
region III by calculating the potential force (—XU) and 
the probability flux (J ss ) (see Fig. [7]). As shown at 
Fig. 0 the potential has local minimum valley along 
the deterministic oscillation trajectory. The values of 
potential are not uniformly distributed along the cycle. 
There is a global minimal potential at a state around 
([p53, [Mdm2]) = ( 1 , 0 . 02 ), which with small curl flux 
forces are large potential forces. The states with larger 
[p53] or [Mdm2] have higher potential, larger curl flux 
forces and smaller potential forces. We note that the 
potential forces are mostly vertical to the cycle path, 
and the curl flux forces are along the cycle. These re¬ 
sults suggest that potential forces tend to attract the 
system states to the oscillation trajectory, while the curl 
flux forces drive the oscillation along the cyclical tra¬ 
jectory. The obvious low potential value at the global 
minimum indicates a pulse dynamics in each cycle the 
system states mostly stay at the state of global mini¬ 
mum while go through the other part of the limit cycle 
quickly. This pulse dynamics is in agree with previous 
studies in mm- 


D. Effects of PDCD5 regulation efficiency 

PDCD5 interacts with the p53 pathway by promot¬ 
ing p53 stability through the disruption of the reaction 
between p53 and Mdm2 nuc . In our model, the p53 degra¬ 
dation rate function is given by a Hill-type function and 
the half-maximal effective concentration (EC50) is de¬ 
pendent on PDCD5 concentration given by Eq. In 
this function, the efficiency of PDCD5 is represented by 
the coefficient aq, larger aq indicates higher efficiency. To 
investigate the effects of PDCD5 regulation efficiency, we 
took (u P 53 , Po) = (0.85, 0 . 8 ) and changed aq from zero to 
triple the default value to examine its effect on system 
dynamics. 

Fig. [ 8^1 shows the bifurcation diagram of p53 with re¬ 
spect to the efficiency aq. There is a supercritical Hopf 
bifurcation (aq = 1.612) as aq increases, where a stable 
limit cycle appears to yield p53 oscillations. The oscilla¬ 
tion amplitude is insensitive with aq. The period increase 
rapidly near the Hopf bifurcation transition, and then de¬ 
creases to 350min when aq increases away the bifurcation 
point (Fig. ^p). 

To further study the stability of p53 oscillation, we 
calculated the barrier height in the potential landscape, 
which is a good measure of the stability of limit cycles 
[33] . The barrier height (BH) is defined by the potential 
function as BH = P max ~ P m in, where U max is the max¬ 
imum potential inside the stable limit cycle, and P m in 
is the minimum potential along the cycle [33] • Higher 
barrier height means more stable. Fig. |§ shows the 
barrier height versus aq. In our simulation, the bar¬ 
rier height is highly fluctuated near the Hopf bifurcation 



%53 


FIG. 5. Typical codimension-one bifurcation diagrams of [p53] with respect to the parameter a P 53 with given PDCD5 level as 
labelled by a-g at Figure [2] In all diagrams, red solid lines represent stable equilibria, black dashed lines are unstable equilibria, 
green solid dots are the maxima and minima of the stable limit cycles, while blue open circles represent the maxima and minima 
of unstable limit cycles. Codimension-one bifurcation point are marked as Fi (or F*) for fold bifurcation points of equilibria, 
Hi for Hopf bifurcation points, LPCi for the fold bifurcation points of limit cycle, HC for homoclinic bifurcation points, and 
SNIC for a saddle-node homoclinic bifurcation point. 


(1.612 < aq < 2.6), and then decreases towards a value 
of about 8.78. 

Next, we considered the effect of aq to the dissipation 
of energy, which measures the inevitability for an open 
nonequilibrium system due to exchanging information 
and energies with its surrounding. The mean rate of en¬ 
ergy dissipation (hd) is give by hd = f F(X) • J(X, t)dX, 
where F is the driving force, and J the probability flux 
vector. According to [42] , the energy dissipation is equal 
to the entropy production rate (e p ) at equilibrium state. 
Therefore, the energy dissipation can be used to obtain 
a global physical characterization of the nonequilibrium 
systems. Fig. [8ji shows the dissipation versus aq. The 
dissipation equals zero when aq is smaller than the Hopf 
bifurcation value and there is a single steady state. This 
shows the detail balance at equilibrium state, in consis¬ 
tent with the fact that the least dissipation cost gives a 
more stable system in general m- When aq increases 


beyond the Hopf bifurcation (aq > 1.612), the energy 
dissipation increases rapidly to maximum at aq = 2.6, 
and then decreases to and sustain at 2 x 10 -7 with aq 
increasing. 


In summary, these results show clearly transition dy¬ 
namics when the system switches from a single stable 
steady state to a stable limit cycle with the increasing 
of PDCD5 efficiency. During the transition, the steady 
state is destabilized to yield p53 oscillations, both oscilla¬ 
tion periods and energy dissipations first increase rapidly 
and then decrease to a sustain level, and the barrier 
height shows large fluctuation in the transition region. 
After the PDCD5 efficiency increases beyond the tran¬ 
sition region, the system dynamics is insensitive to the 
efficiency. 
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FIG. 6. The potential landscapes for the ten typical dynamics. Parameters are the same as in phase diagrams of Fig. [3] and 
the diffusion matrix was taken as the diagonal matrix with noise strength D — 1.0 X 10 -6 . Refer to the text for details. 


IV. CONCLUSIONS 

The dynamics of tumor suppressor p53 plays important 
roles in the regulation of cell fate decision in response to 
various stresses[44]. PDCD5 is known to interact with 
the p53 pathway and functions as a co-activator in p53 
regulations. In p53 pathways, the p53-Mdm2 oscillator is 
crucial for cell response to DNA damage. In this study, 
we have systematically studied effects of PDCD5 to the 
p53-Mdm2 oscillator by methods of bifurcation analy¬ 
ses and potential landscapes. Our results reveal that the 
p53-Mdm2 oscillator can display monostability and bista¬ 
bility under low PDCD5 expression. When PDCD5 level 
is upregulated (for example, by DNA damage), p53 oscil¬ 
lations emerge by Hopf bifurcation, and bistability with 
the coexistence of both stable oscillation and a stable 
steady state is possible for proper PDCD5 level and p53 
synthesis rate. These results were further verified by the 
potential landscapes, which clearly show the transition 
of landscapes with changing parameter values. We have 


also investigate the effects of PDCD5 efficiency in the 
interaction with p53 pathway, we showed that p53 os¬ 
cillations can only be induced only when the efficiency 
is larger than a critical values of Hopf bifurcation, and 
the system dynamics show clear transition features in 
both barrier height and energy dissipation when the ef¬ 
ficiency is close to the bifurcation point. Such abnormal 
behaviours at the transition region have been highlighted 
in recent years in the application of prediction and early 
diagnosis of complex diseases [45H47]. 


In this study we have focused on the PDCD5-regulated 
p53 dynamics in the p53-Mdm2 oscillator. Our results re¬ 
veal global p53 dynamics mediated by PDCD5 and the 
levels of p53 production. However, a more complete p53 
pathway is certainly necessary in understanding cell fate 
decisions in response to DNA damage. Further consider¬ 
ation of the effects of PDCD5 on a complete p53 network 
is certainly required in future studies. 
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Hopf bifurcation point, red solid lines represent stable steady 
states, black dashed lines are unstable steady states, green 
solid dots are the maxima and minima of stable limit cycles, 
(b) Period of stable limit cycles, (c) The barrier height of sta¬ 
ble limit cycles (BH). (d) The energy dissipation (e p ). Here 
u p 53 = 0.85, Pq m 0.8, and the noise strength D — 1.0 x 10 -6 . 
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